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Abstract 

In this paper we construct the noncommutative Grosse-Wulkenhaar model on 2- 
dimensional Moyal plane with the method of loop vertex expansion. We treat renor- 
malization with this new tool, adapt Nelson's argument and prove Borel summability 
of the perturbation series. This is the first non-commutative quantum field theory 
model to be built in a non-perturbative sense. 
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1 Introduction 

Quantum field theories on noncommutative space time became popular after the discovery 
that they may arise as effective regimes of string theory either due to the compactification 
[1] or due to the presence of the Green-Schwarz B field for open strings [21 [3] . The simplest 
non-commutative space is the Moyal space, which could be considered also as a low energy 
limit of open string theory. However the quantum field theories on the Moyal space are 
non-renormalizable due to a phenomenon called UV/IR mixing, namely when we integrate 
out the high scale fields for non planar graphs, there exist still infrared divergences and the 
divergent terms cannot be compensated by counterterms 

Several years ago H. Grosse and R. Wulkenhaar made a breakthrough by introducing a 
harmonic oscillator term in the propagator so that the theory fully obeys a new symmetry 
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called the Langmann-Szabo duality. They have proved in a series of papers O El [7] (see 
also p]) that the noncommutative field theory possessing the Langmann-Szabo duality 
(which we call the GW4^ model hereafter) is perturbatively renormalizable to all orders. 

After their work many other QFT models on Moyal space [HI HDl IHl HH 1121 1131 US] or 
degenerate Moyal space [161 HZ] have also been shown to be perturbatively renormalizable. 
More details could be found in [20l [211 122] . 

The GW4^ model is not only perturbatively renormalizable but also asymptotically safe 
due to the vanishing of the (3 function in the ultraviolet regime [21], [2S], [2S], and the 
renormalization fiow of the coupling constant is bounded, which means that it is even better 
behaved than its commutative counterpart. In the commutative case we either have the Lan- 
dau ghost for the (pf theory and QED in the ultraviolet regime or we have confinement for 
non-Abelian gauge theory in the infrared regime. This makes the non-perturbative construc- 
tion of these commutative field theories either difficult or impossible; it also means that GW4^ 
is a prime natural candidate for a full construction of a four dimensional just renormalizable 
quantum field theory |28j . 

Constructive field theory or constructive renormalization theory builds the exact Green's 
functions whose Taylor expansions correspond to perturbative quantum field theory. The 
traditional techniques for Bosonic constructive theories are the cluster and Mayer expansions 
which divide the space into cubes and require locality of the interaction. However due to the 
noncommutativity of the coordinates in Moyal plane, the naive division of space into cubes 
seems unsuited in this case. What's more, due to the non-locality of the interaction vertex, 
it is not clear to which cube an interaction vertex should belong. So it seems very difficult, 
if not impossible, to construct the Grosse-Wulkenhaar model with the traditional methods 
of cluster and Mayer expansion. 

In this paper we shall construct the Grosse-Wulkenhaar model in the 2 dimensional 
Moyal plane or GW2, as a warm-up towards building non-perturbatively the full GW4 model. 
The method we use is the Loop Vertex Expansion (LVE) which was invented precisely for 
overcoming these difficulties of the traditional methods [27] and has been applied also to other 
QFT models with cutoffs [^ EHl EI] • This method has been tested first on the ordinary </)| 
commutative case in a companion paper [32] ■ Another approach towards the construction of 
GW4, based on the combination of a Ward identity and Schwinger-Dyson equations, is given 
in [33]. A non-perturbative numerical study of the noncommutative quantum field theory 
model is given in [18] . 

The loop vertex expansion is a combination of the intermediate field techniques with the 
BKAR forest formula [321 SD] • In the case of a interaction, the first step is to introduce 
an intermediate field a and integrate out the original field 0. After that we get a non local 
interaction (l/2)Tr log(l + zyXC^/VCVS)^ where Tr means summing over all indices for the 
matrix model (or integrating over the whole space for the commutative case) while keeping 
the cyclic order, and C is the covariance of the original fields (p. After this step we have a 
new representation of the same theory where the dynamical variables are the intermediate 
fields a. The new features are that the interaction vertex becomes non-local (even when 
the initial vertex, like in ordinary 0^ was local) while the covariance of the a fields becomes 
ultralocal, which means that the a fields have trivial propagation. 

Then we apply the BKAR tree formula to obtain the (unrenormalized) connected Green's 
functions. But this is not the whole story. Some parts of the Wick ordering counterterms 
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can be easily compensated with the corresponding terms within the LVE formahsm. But the 
Green's functions obtained after the LVE are still divergent due to other remaining counter- 
terms from the Wick ordering of the interaction which correspond to terms which are not 
explicitly visible within the LVE. This is why we need a second expansion to generate terms 
called inner tadpoles and to compensate them with their associated counter-terms. We 
call this second expansion the cleaning expansion; it was introduced and applied to the 
commutative 02 case in [32] . 

After this cleaning expansion, using analogs of the Nelson's bounds for 0| (see [21] and 
references therein), we obtain convergence of the expansion hence rigorous construction of 
the Green's functions and we can prove their Borel summability. In fact due to the harmonic 
Grosse-Wulkenhaar potential these bounds are easier in the GW2 case. 

This paper is organized as follows. In section 2 we introduce the basic properties of the 
noncommutative Moyal space, the matrix basis for functions defined on this space, the GW2 
model in the matrix basis and the combinatorial factors for its Wick ordering. In section 

3 we introduce the intermediate field techniques and the BKAR tree formula. In section 

4 we consider the amplitude, the graph representation and the dual representation of the 
connected functions obtained by the LVE. In section 5 we define the cleaning expansion. In 
section 6 we consider the renormalization. Sections 7 and 8 are devoted to the resummation 
of the counter-terms and the Borel-sumability. Two examples for the first and second order 
perturbative renormalizations are given in the Appendix. 



2 Moyal space and Grosse-Wulkenhaar Model 
2.1 The Moyal space 

The D-dimensional Moyal space for D even is generated by the non-commutative co- 
ordinates x'^ that obey the commutation relation = iO^^", where Q is a D x D 
non-degenerate skew-symmetric matrix such that G'^'' = — G'^'^. It is the simplest and best 
studied model of non-commutative space (see [121120] for more details). 

Consider the Moyal algebra of smooth and rapidly decreasing functions iS(]R^) defined 
on M.^. S{M.^) is equipped with the non-commutative Groenwald- Moyal product defined by: 
yf,geSD:=SiR^), 

1 



{f*eg){x) = I ^-^d^yf{x+\e-k)g{x + y)e'^-y 



TT^IdetGI 



d''yd''zf{x + y)g{x + z)e-''y''-'^ . (2) 



We now set D = 2 in this paper. We could define the creation and annihilation operators 
in terms of the coordinates of the Moyal plane [6], IT] [19] : 

a = -j={xi + 1x2) , a = —^{Xi - 1x2) , 

V 2 V 2 
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such that for any function / G Sd we have 



01 df 6i df 

(a^/)(x) = a{x)f{x) + -jg^i^) , = a{x)f{x) - -^Q^i^) , 

(a^/)(x) = -a{x)f{x) - ||^(a;) , (/^a)(x) = a(x)/(x) + ||^(^) . (4) 

With the creaton and annihilation operators we could build the matrix basis for the 
functions 0(x) G iS(]R^) defined on Rq, such that: 

^(x) = ^(Pmnfmn{x), (5) 
mn 

where the matrix basis fmn{x) is defined as: 

fUx) := I , a^" ^ /o ^ a^'^ (6) 

mm(m,n) / \ / \ 

and /o(x) = 2e"5(^i+^2). 

The matrix basis defined above has the following properties: 

(/mn * fkl){x) = 6nkfml{x) . (7) 

and 



2.2 Ribbon graphs 

Matrix models or quasi-matrix models use ribbon graphs, which were first introduced by 't 
Hooft for QCD [36j. A ribbon graph G = {V,E) is an orientable surface with boundary 
represented as a set of closed disks, also called vertices, V and a set of ribbons, also called 
edges E, such that all disks are attached to each other by the ribbons glued to their bound- 
aries. Figure [T] is an example of a non-planar ribbon graph. The interested reader could go 
to EH] for a more general introduction. 




Figure 1: An example of a ribbon graph. 
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2.3 The 2-dimensional Grosse-Wulkenhaar Model 

The 2-dimensional Grosse-Wulkenhaar model or GW2 model is defined by the Wick-ordered 
action: 

+ (9) 

where (j){x) is a real function defined over the noncommutative coordinates x^, = 2(0^^)^,^x'' 
and the Euclidean signature has been used. Remark that this model is not translation in- 
variant due to the harmonic oscillator term. 

The action ^ has a remarkable symmetry such that if we exchange the position and 
momentum operator: 

^ Xf,, (f){p) ^ 7rv^|det^|0(x), (10) 

where 

0(p) := J £xe-'P®''^{x), pBx ■.= p^e>"'x^, (11) 
the action ^ is invariant up to a scalar factor: 

S[<P;f,,X,Q]^^'S[<P;^,^,^]. (12) 

This symmetry is called the Langman-Szabo duality [23]. It is essential for solving the UV/IR 
mixing problem hence for proving the renormalisabiliy of the Grosse-Wulkenhaar model. 

The renormalisability of the GW model has been studied extensively in the matrix repre- 
sentation [71 IB] , the direct space representation [3l] and the parametric representation [35] . 
In this paper we shall use the matrix representation. 

In the matrix basis the GW2 action could be written as: 

^[0]=27r^Tr[^0A0+:^0l 
1 , . , A 



2tc9 ^ ^-(pmn^mn.kKt^kl + ^ : <Pmn(t)nk<Pkl<Plm ■ ■ (13) 



m,n,k,l 

where (pmn is a Hermitian matrix. In this basis the Laplacian reads: 

2 2 

^mn,kl = [/U^ + ^(m + n + l)]Sml5nk " ^(1 - (14) 

X [a/ (m +l){n + l)6m+i,iSn+i,k + Vmn6m^i^iSn^i,k] ■ (15) 

We introduce also an ultraviolet cutoff A for the matrix indices such that < m, n < A, 
then the covariance matrix Cmn,ki is a A^ x A^ dimensional matrix in general. We could 
obtain the covariance Csr,ki by the following relation [6|: 

^ ^ ^mn,rsC sr,kl ^ml^nk- (■^^) 

r,seA 



When ^1 = 1, the kinetic matrix reduces to a much simpler form: 

r 4 1 

^mn,kl= fJ^"^ + -^{m+n+1) 6mlSnk, (17) 

which is a diagonal matrix in the double indices and the covariance matrix Cmn,ki reads: 

1 

Cmn,kl = ; J7 ; "TT'^m/'^nfc- (l^) 

Since = 1 is a fixed point of this theory ^25j, we shall for simplicity take f2 = 1 in the 
rest of this paper and write the covariance as Cmn = Cmn,nm = m^n+i ^'^^ simplicity. 
The Wick ordering of the interaction with ultraviolet cutoff A reads: 



>mn(Pnk(Pkl(Plm ■— <Pm.n<Pnk<Pkl<Plm — 0<Pmp<Ppm 



T^ + 6Tr^(r^)^. (19) 



where 



and 



m 



= log ~ log A, for 1 ^ m < A, (20) 

n 4- m m 



q + m m 

9=0 ^ 



n = Tr^(r^)^ = ^)(E ^) - (21^'2 + ^)A. (21) 



m p q 



The counter terms Tr(/)^T and and an explanation of the combinatorial factors are 
shown in Figures |2] and [3j 




Figure 2: The counter-term Tr(0^T) from the Wick ordering. There are 4 possibilities of 
contracting the nearest neighbor fields, for example, 1 with 2, 2 with 3 , 3 with 4 or 4 with 1; 
And each contraction has additional two possibilities of whether the two uncontrcted fields 
are in the external face of the inner face, as shown in this Figure. In total this gives a factor 
8. 

Note that since the infinite volume limit is hidden into the matrix basis, we don't need 
any more the cluster expansion. In this sense the Grosse-Wulkenhaar model is even simpler 
than the commutative (f)\ theory. 
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2 



3 




Figure 3: The counter term and combinatorics. There are 2 possibihties of forming a 
vacuum graph, for example 1 contracts with 2 and 3 contracts with 4, or 1 contracts with 
4 and 2 contracts with 3. For each case there are an additional 3 possibilities shown in this 
Figure. In total this gives a factor 6. 

3 The intermediate field representation and the Loop 
vertex expansion 

3.1 The intermediate field representation 

The partition function for the matrix model reads: 

Z{X) = [ d^i{<Pmn)e-'^^-\ (22) 



where 

rf/i(0™„) = 7r-^'/2e-i/2E™<^™^™^-] n^Re(</'™-)^Im('/'--). (23) 



is the normalized Gaussian measure of the matrix fields (j)mn with covariance C given by (18) 
and S[(l)m-n\ is the Wick ordered interaction term. 

We introduce the Hermitian matrix intermediate field and the partition function 

can be written as: 



Z{X) = J d^{a)d^{(j))e~^^'^'^*''-'^'^"^'^'^''''-'^'''-~^'^^^*'^^ (24) 
where 

f//i(a„„) = 7r-^'/2g-i/2^™'^- J]rfRe(a^„)Im(a„„). (25) 

mn 

is the normalized Gaussian measure for the hermitian matrices a with the covariance: 

= i„A.,. (26) 

Integrating out the matrix fields 0mn we have: 

Z{X) = J c//i(a)e^*^'^'^'''^"^T'^°s[^+*^ ci/2(/®a+a®7)ci/2 j+fAT^^ ^27) 



where TrT^a = X^m-^m) ^ is the A x A dimensional identity matrix. It is useful to keep all 
the indices of the covariance C = {Cmn,ki} and write is as a A^ x A^ dimensional diagonal 
matrix with elements Cmn,ki = m+n+i ^"^' ® ^'^k- C^^"^ is the square root of the covariance C 
with elements: 



[C^^'^]mn,kl — \ ; —r^ml ® ^nk- (28) 

m + n + 1 



The term Trlog[l + iy2\ C^^'^{I ® cr + a ® I)C^^'^ ] is called the loop vertex, which was 
first introduced in [27]. Remark that the term ® cr + o" ® /]C^/^ in the loop vertex is 

Hermitian. Graphically this term means that the intermediate field a can attach to either 
the inner border or the outer border of a ribbon graph. 

We define the A^ x A^ dimensional matrix a as 

a = I(g)a + a(g)I, (29) 



then the loop vertex reads: Trlog[l + i^/2X C^I'^dC^I'^ ]. In the rest of this paper we shall 
write (J as (T for simplicity and the loop vertex reads: Tr log[l + i\/2\ C^I'^oC^I'^ ]. We just 
need to remember that C^l'^ and a are A^ x A^ dimensional matrix. 

Before using the BKAR tree formula to calculate the connected function, we first of all 
expand the loop vertex term as: 

-^Trlog[l + zy2AC^/VC^/2] 

= -^Tr(\/2AC^/VC^/2) - ^Trlog2[l + %^f^C^I''oC^I\ (30) 

where log„(a;) for n > 2 is defined as the n-th Taylor remainder of the function log(x): 

log„(a;) = log(x) - [x - 12 + j'i ■ ■ ■ + (-l)"+^x"/n]. (31) 



The first term in formula (30) reads: 

A 

2 



jTr(v^Ci/VCi/2) = -2 ■ -JYxY^T^^amm. (32) 



where the additional factor 2 is due to that there are two kinds of tadpoles (see Figure |4|) 
when taking the Trace after using the fusion rule for the covariances 

Ci/2. Hence the leaf 



Figure 4: Tadpoles. 



tadpole term (32) could cancel partially the linear counter term of the Wick ordering, see 
formula (27). 
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So we could define the new interaction vertex as 

V{a) = -^TT\og2[l + tV2XC^/^aC^/^] 

+ tV2XTTT^a+^XTl (33) 

where the leaf tadpoles have been already canceled by the counter term and the partition 
function could be written as: 

Z{X) = J di2{a)e^^''\ (34) 
3.2 The BKAR Tree formula and the expansion 

The most interesting quantities in our model are the connected Schwinger's functions. We 
shall derive them by the BKAR tree formula. 

Let us first of all expand the exponential as ^^^-^ . To compute the connected func- 
tion while avoiding an additional factor n\, we give a kind of fictitious index v, v = 1, ■ ■ ■ ,n 
to each field a in vertex V{a) and we could rewrite the expanded interaction term as 
X^n n"=i ^^n\ ^ ■ This means that we consider n different copies (T„ of a with a degener- 
ate Gaussian measure 

n 

di^iWv}) = duicTvo) Yl ^^^'v ~ '^vo)da^' , (35) 

where Vq is an arbitrarily marked vertex. 

The vacuum Schwinger's function is given by: 

Theorem 3.1 (Loop Vertex Expansion [27]). 

oo ^ 

logZ(A,A,V) = ^- (36) 

n=l T with n vertices 

= { n E [ /' dwi] I / duriia^}, {w}) 

i£T 0(Tmt,ni (J'^pt,qe J v=l 



where 



each line i of the tree joins two different loop vertices V'"^^'^ and which are iden- 

tified through the function Smig'^Sn^p'^,, since the propagator of a is ultra-local. 

the sum is over spanning trees joining all n loop vertices. These trees have therefore 
n — 1 lines, corresponding to a propagators. 
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the normalized Gaussian measure (iz/7-({cr^}, {w}) over the fields a^, has now a covari- 



ance 



< ^mn. ^kl >= SnkSmlw'^iv, v' , {w}) , (37) 

which depend on the "fictitious" indices. Here w^{v,v', {w}) equals to 1 if v = v' , and 
equals to the infimum of the We for i running over the unique path from v to v' in T 
if V ^ v' . 

The proof is a standard consequence of the BKAR formula (see [271 129]). 

The n-point Euchdean Green's functions or the Schwinger's functions, with n = 2p even, 
could be derived similarly by introducing the (matrix) source fields j = [jmn] ■ The partition 
function with source fields reads: 

Z{A,J{X)) = j d;,(0)e-lTr[(0^-4TA)2+lOTl]+Tri<^ (gg) 
^^(^^)g|AT2+Tr(2iv^TAa-|log(l+*V2ACi/2aC'i/2)) + lTr[jCi/2/?C'i/2j] 



where 



is called the resolvent matrix. We define also the full resolvent matrix: 



R = = (39) 

1 + iy2ACVVCV2 



Ria) = C'/' =^ ^1/2. (40) 

1 + zV2ACi/2aCi/2 



The connected n-point function is given by: 

1 d'^P 

S{hru--- , V2p) = 7^7^^ logZ(A,j)li=o. (41) 

More precisely, we have the following theorem for the connected Schwinger's function: 
Theorem 3.2. 



TT n=0 T with n vertices ^ leT ^ ^ 



n 2p 

n n^W). (^)' (42) 

v=l k=l 

where the sum over tt runs over the parings of the 2p external variables according to the 
cyclic order and -R«i^(fe),ri ^(j., (o") o'^e the resolvent matrices whose explicit form are given by 



formula (39). 



So in this paper we consider only the vacuum connected Schwinger's function, as the 
n-point function could be easily obtained by derivation w.r.t. the source fields. 
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4 The graph representations and the amplitudes of the 
LVE 



4.1 Direct representation of the LVE 

Now we consider the basic graph structure of the Loop vertex expansion, which we call the 
direct graphs or the LVE in the direct representation. 

There are three basic line structures or propagators in the LVE: 



The full resolvent R{X,a) (see also (40)) is defined as: 



R can be written in terms of the resolvent matrix defined in (39) as: 



R^n = [C'/' R C'^Xn- (44) 



• The propagators Cmn between the original fields (^mn, 

• The propagators between the a fields. 
The propagators are shown in Figure |5| 



Figure 5: The propagators in LVE. A stands for the resolvent Rmn, B stands for the pure 
propagator Cmn, and C is the propagator of the a fields. 

There are three kinds of interaction vertices in the LVE: the counter terms, the leaf terms 
K with coordination number 1 and the general interaction vertices V. A leaf vertex K is 
generated by deriving once w.r.t the a field on the logg term: 



-lV2\\C^'\ =^ l)C^''']ran. (45) 



2 l+iV2ACVVCV 

A general loop vertex could be obtained by deriving twice or more with respect to the o 
fields: 

K^.„,_,r.,(A,a) = — ^ [-^Trlog(l + zv/2ACi/VCi/2)] 

= -^-(tV2Xn-ir J2 R^.m^,., A) ■ ■ ■ Rm^^^.^m. (a. A). (46) 

r 
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A 



B 



C 



Figure 6: The basic graph elements of LVE. Graph A means the leaf K, graph B means the 
counter term and graph C means the most general loop vertex which has several a fields 
attached. 

with p >2 and the sum over r is over the p cyclic permutations of the resolvents. 



The basic interaction vertices are shown in Figure [33} where we didn't show explicitly 
the pure propagator C^/^. 

After integrating out the a fields we obtain the vacuum graphs which are planar ribbon 
graphs without broken faces. The first order graphs whose amplitudes are proportional to A 
are shown in Figure [Tj Graph C in this figure means the counter term |AT^ that appears 
only in the first order expansion. 




Figure 7: The graphs of first order expansion. We have taken account of all possible config- 
urations of graph A in the calculation. 




Figure 8: A general graph of Loop vertex expansion. Here the leaf means the leaf loop vertex 
K. 

Since we consider a real symmetric matrix model, the leaf loop vertices and the counter 
terms could appear either in the outer side, in the inner side or both sides of a loop vertex, 
we call the graph of the first case the normal direct tree and the last two case the inner direct 
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tree. The normal direct tree and the inner direct tree are shown in Figure |8] and Figure [9} 
respectively. 

The amplitude of the two graphs are exactly the same, so in this paper we only consider 
explicitly the first case for the direct graph representations and the dual representations. 
The contributions of the inner direct trees to the amplitude should be taken into account 
by calculating correctly the corresponding combinatorial factors. A general graph of LVE is 




Figure 9: A tree of loop vertices where the counter terms or leaf terms are attached to the 
inner part of the loop vertices. Here the leaf means the leaf loop vertex K. 

shown in Figure |8] and Figure |9] 

Now we have the following definitions for the objects in T: To each counterterm TTaT^D^r 
where Dmn is the pure propagator generated by the derivation on the loop vertex, we 
associate a tree line index ic] to each leaf term Kmn{o') we associate a tree line index i and 
to each full resolvent R{a)mn sandwiched by the two tree lines ii and ir, such that ii is on 
the left side and ir is on the right side, we associate two tree line indices £i and ir ■ 

So the amplitude for a tree of loop vertices reads: 




where the global Tr means that we need to take the cyclic order of all the objects. The 
indices £; or ir may coincide with the indices i or ic- 

4.2 The dual representation 

Since a LVE graph in the direct representation is planar, the notion of duality is globally 
well-defined. In the dual representation we have a canonical (up to an orientation choice) 
and more explicit cyclic ordering of all ingredients occurring in the expansion [32j . 

^ Remark that for the orientable model with interaction X4>(j>(f>(p the loop vertex has the form Trlog(C(TL(8) 
1 + 1 (g) 1) or Trlog(l (K) 1 + 1 (E" <JrC), which means that the a fields appear only on one side of the loop 
vertex (left or right). See also [27] for a more detailed discussion on this point . 
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Figure 10: The dual graph of a LVE. The area enclosed by the bold dash ribbons correspond 
to the original loop vertices. 



. In the dual representation the original tree of loop vertices T is replaced by the dual 
tree T, which consists of a single huge loop vertex C with dual tree line I. Each dual tree 
line £ G T corresponds to a line £ G T in the direct picture, and is pictured as a bold dash 
line. See Figure [TO] for a dual tree. 

The important fact reflecting the tree character of T is that these bold dash lines cannot 
cross. Hence they divide this disk C into different connected regions v &V , each of which 
corresponds to a single loop vertex f G of the direct expansion. This explains how the 
two pictures, dual of each other, are equivalent and how anyone can be reconstructed from 
the other. 

This cycle C contains all the ordinary full lines (0 resolvents) of the direct representation 
and the counter terms, but reread in the cyclic order obtained by turning around the tree. 
The order of all the objects in the dual representation should be the same as in the direct 
representation. 

Here is the notation for the objects (9 = i?US'UTu£in the dual representation: 

• The counter-terms that are pictured as black vertices b E B. More exactly, each 
counter-term is considered as a black vertex attached to a pure propagator, so that 
the amplitude of each counter-term equals to the product of the amplitudes of the two 



objects. See formula (48). To each counterterm we associate an index 



The bold dash lines £ G T that correspond to the original a propagators. We also call 
them the dual tree line. 

Simple leaves Vq G S, i.e resolvent hues surrounded by a single dual tree line £. For 



example, the normal lines in region A and E in Figure 10 correspond to two simple 
leaves. To each leaf we associate an index i. 

Non-leaf resolvent lines Vi = {ii,ir} G £ , ie resolvent lines not in 5* that correspond 
to the full resolvents R^^''{a) in the direct representation. We follow the clockwise 
direction in the dual graph and use two indices ii and Ir to indicate each non-leaf 
resolvent line sandwiched between the two different dual tree lines ii and ir- 
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The beauty of the LVE representation is that all the various traces of the loop vertices 
in the direct representation give rise to a single trace in the dual representation. This is 
the fundamental observation which made the representation suited for constructive matrix 
madel and non commutative quantum field theories [27j. 

We need a label u to describe the various objects met when turning around the cycle C. 
Then to each object u E O is associated an operator Pu, with value 

P„ = Ta^L if u = beB, (48) 
Pu = Ki^a) a u = v^eS, (49) 
Pu = B^±{(y) a u = vieC, (50) 

In this dual picture the measure dv{a, w) corresponds to the following rule: the weakening 
factor between a and a a^)/ is the infimum of the w parameters of the lines i that have to 
be crossed to join the two regions v and v'. 

Hence the amplitude for a tree of loop vertices reads: 

PA = logZ(A,A,V)= ^ Gr 

T lines, |T|=n— 1 

Gr={-Xr^'l rfK^,^)Tr|n^^^P„}. (51) 

where the operator Tr means that we sum over all matrix indices following the cyclic ordering 
of all the objects P„. 

Theorem 4.1. The vacuum Green's function Pa is absolutely convergent and defines an 
analytic function in the half-disk "Da = I 3ft(A^^) > log A, where K is a large constant. 

Proof. This is just an application of the techniques of [271 129] which we recall briefly. 

• The number of trees is n"~^, by Cayley's theorem, 

• each full resolvent R{a)mn is bounded by 1 in the disk, 

• each decorating counter term is bounded by log A. 
So the connected function is bounded by 

Pa < y ^A"(lnA)" < V (Air)"(ln A)". (52) 
^-^ nl ^ — ^ 



n 



which is a convergent geometric series as long as 3ft(A~^) > i^'logA but divergent when 
A — oo. So the loop vertex expansion alone couldn't tame all divergences and we still need 
renormalization. Remark that the role of the divergent tadpole terms and counter-terms are 
in some sense changed. The divergences are coming from the linear counter terms due to the 
Wick ordering of the interaction and should be compensated by the inner tadpoles, crossings 
and nesting lines (see page 19 and Figures (13), (16) for the definitions and pictures of these 
objects) that we generate by contracting the a fields hidden in the resolvents. 

So that we shall put the perturbation series and the remainders on the same footing, 
which is another interesting property of the loop vertex expansion. We shall discuss this 
point in more detail in the next section. □ 
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5 The Cleaning Expansion 



In this section we shall introduce another expansion, which we call the cleaning expansion, 
to compensate the divergent counter terms. We shall work in the dual representation of 
the LVE, since in this representation the cyclic ordering of different objects is more explicit. 
Instead of having a trace operator for every objects, we have a single global trace operator. 
We shall use the multi-scale representation of the propagators and the resolvents. 

We impose also the stopping rule for the cleaning expansion to make sure that we don't 
expand forever (this would lead to divergence). We stop the expansion until we have gained 
enough convergent factors to compensate the divergent Nelson's factor. 

5.1 The sliced propagator and resolvents 

We introduce the Schwinger parameter representation the propagator: 

1"^ 2 /"I 2 

C^n = / rfae-"('^+™+"+i) = K rfae-°('^+"^+"+^^ (53) 
Jo Jo 

We decompose the propagator as: 

oo 

Cmn = ^ C'mn) (^4) 
3=0 

,^-(2i-2) 

Ci^ = / ciae-"('^'+^(™+"+^)\ (55) 

and M is an arbitrary positive constant. We could easily find that 

\CU < i^M-2^e-*^"'ll^'+t(™+-+i)ll. (56) 

We have the associated decomposition for the amplitude of the connected function: 

Gr = J2G'T- (57) 

j 

Due to the cyclic order of the global Trace operator, we could rewrite the loop vertex 



where 



in the non-symmetric from Trlog[l -|- i\/2XaC] and the resolvent defined in formula (39) is 
written as: 

Rmn = -7, [--Trlog(l + 2V^crC)] = V[ ^^^]„pCp„. (58) 

OCTmn ^ ^ ^ l + ZV2AcrC 

It is convenient to introduce further notations to write down in more compact form the 
contribution, or amplitude of a tree in the LVE. We put 



-Dmn — W 2ACmn; (59) 
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so we have 

R{a) = — and R{a).mn = [R{a)DUn, (60) 

i + crIJ 

where Dmn and R{(j) should be thought of as a matrix operator. 
Then we define the shced resolvent as: 

R'i'^) = = T^, (61) 

and the full resolvent as 

Rl^ = [D^ R^Un. (62) 
For each resolvent we have the algebraic equation 

R^{a) ^ ^ 



jk<j ^ Z-^k<j 

= ^ p,-l , = - R'-\o)aDm^. (63) 

1 + R^ '-[ajaD^ 

5.2 The cleaning expansion 

In this section we shall expand the resolvents and contract the a fields to generate the inner 
tadpoles for canceling (partly) the counter terms and to generate the crossings and nesting 
lines for gaining convergent factors. This is called the cleaning expansion [52] . 

The canonical cyclic ordering of the LVE is also essential for the cleaning expansion, as 
otherwise it would be difficult to find the corresponding "weakening parameters" rule for 
that Taylor remainder term. But fortunately the cyclic ordering solves nicely this problem! 
The beginning of the cycle C is explicitly and fully cleaned of potential tadpoles; the rest of 
the cycle C has no weakening parameters on the remaining potential tadpoles! 

We won't just expand each resolvent naively as power series since that 

would generate unnecessarily complicated weakening factors when we contract the a fields of 
different regions in dual graph, and what's more, this naive expansion is not convenient for 
the multi-scale analysis of the amplitude, which is crucial to analyze the Taylor remainder 
terms. 

We apply the cleaning expansion in two steps: 

In step 1, we start from an arbitrary marked point in the dual graph and apply the 
algebraic induction formula ([64j) 



Rj(a) = R^-\(r) - R^-\(t) a (64) 

which rewrites the multi-scale representation of each resolvent of scale j into two parts, for 
the resolvent we meet clockwisely. The first part is resolvent of a lower scale j — 1, and the 
second part is a product form of a pure propagator of scale j, a resolvent of scale j and a 
resolvent of scale j — 1. We start form the highest scale jmax- Due to the cyclic ordering the 
pure propagator term appears only on the left side of the product. 



So by using the formula (64) inductively we could lower down the scale of the resolvents 
and propagators on the one hand and decompose the resolvents into a cleaning part which is 
the pure propagator with a field attached and the uncleaned resolvents, on the other hand. 
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In step 2, we use the integration by parts of the a fields with respect to the Gaussian 
measure. For example, for each term in the second part of formula (64) we have: 



y d/^(cr)i?-^""'" ^(o-)o-L)-'i?-'™""(o-) X other resolvents of scale jmax 
r d ■ 

= dn{a)- [i?^'"''""^((T)D-''"""i?^™"((T) X other resolvents]. (65) 

J OCnm 

We forget the case where ^ acts on the .R-J™"^-^ term as this gives minor contribution and 
we consider only the case ^i?-^™"^ for the first resolvent i?-?™-"^ we meet, for which we have: 

= -D^ W^D^]R^. (66) 



Here we have ignored the matrix indeces. This corresponds to the second line of Figure(ll), 
which shows the typical cleaning expansion process. Then for each resolvent sandwiched by 
two pure propagator D we use the formula = 1 — x^^: 

j 

= Smn-J2iJ2''^'^-^P^in- (67) 



The first term in formula (67) means we identify the two pure propagators generated from 
the cleaning expansion so that we generate an inner tadpole. Remark that from the above 
process we don't generate a propagator with a single scale but ^j™^.^ , the tadpole 
we generated means the sum T\ = ~ In A. This corresponds to the left hand side 

of the third line in Figure (11). In the second term of formula (67) we have again a single 
a field in the numerator and we shall perform the integration by parts. There exists still 
several possibilities: 



• If the derivation ^ act on the nearest resolvent, we shall repeat the process (67). 

• if the derivation act on the second resolvent, this would generate a crossing line. 

• it is possible that the a always contracts with the nearest resolvents without generating 
a crossing nor an innertadpole. In this case the lines of a propagators are called nesting 
lines. See Figure (12) . 

We ignore the counter terms which could be anywhere in the dual graph, just for better 
explanation, as they don't change the cleaning process. 

Now we consider the amplitudes for crossing lines and nesting lines. A crossing means 
a crossing of two cr propagators that correspond to 6 functions and are pictured as dash 
lines, and the a propagator that overlaps another is attached to a pure propagator. Since a 
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Rj Rj"' Rj"' ioj Rj 

Rj"' ;'dJ Rj \Dj R-' 



I ' ' I - + I 'I 'I 

dJ rJ 

Figure 11: The cleaning expansion. The ordinary ribbon stands for the resolvents and the 
bold double hne stands for the pure propagator. The dashed lines should be envisioned also 
as double hues and these stand for the a fields, u is the marked point. 

crossing line doesn't form a tadpole, the amplitude of each crossing line equals to the bound 
of the sliced propagator (see formula ([56|)): 

IIC^^II ^ |irM-2^e-*^"''ll^'+5('"+"+^)ll| < KM-''\ (68) 

Nesting lines means many a propagators that each one is overlapped by another. Each 
nesting line is made up of a a propagator with two pure propagator attached to each end 
of the a propagator, and many (at least one) resolvents or /and pure propagators being 
overlapped by the a propagator. 

The amplitude of a nesting line is bounded by: 

l^if iH^S'll < l^it'CU ^ < KM-'^, (69) 

where [Y[ R^m^n''^] means the full resolvents overlapped by the nesting line that we consider 
and means the pure propagator attached to either side of the a propagator. Since they 

don't form any closed tadpole, it doesn't matter which C4„ we consider; we have used the 

i' I' j 

fact that each full resolvent R^i^' is bounded by 1. 

So we gain a convergent factor M"^-' for each of the two cases. 



Figure 12: A graph with 2 nesting lines. 
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u 




Figure 13: The crossing generated by the case that the field cti is contracted with other 
resolvents which is not the first one in formula (65), for a more general case; The bold 
double dash line means the original dual a propagator, the normal double dash line means 
the a propagator generated by the cleaning expansion. Clearly this forms three crossings, u 
is the marked point. 



5.3 The stopping rules 

We cannot do the cleaning expansion forever, as this would develop the full perturbation 
series and ultimately diverge. More exactly, for each scale j the number of crossings shouldn't 
exceeds n(j) = (aj)!, where a is a positive constant to be fixed, otherwise we would no longer 
be able to compensate for that divergent factorial with the convergent factors coming from 
the crossings and/or nestings. 

So we impose the following stopping rules for the cleaning expansion and write the rest 
terms as Taylor integral remainder when enough cleaning, depending on the scale, has suc- 
ceeded: 

• We start from resolvents and pure propagators of scale jmax and go on the expansion 
towards lower scales. 

• We stop the expansion whenever we generate an inner tadpole and we cancel it with 
its corresponding counter-term. 

• We go on the expansion when we generate a crossing or nesting lines until the total 
number of crossings and nestling lines reaches n{jmax) = CLjmax B ^^^^ total we 
could gain a convergent factor M""-^™"^ from the crossings and the nestling lines while 
lose a divergent combinatorial factor (ajmax)! ~ e^-^™"^^"-'"''"', which could be easily 
compensated. 

• It is very possible that we couldn't generate enough crossings by expanding a single 
resolvent in the dual tree. Then we need to go to the next resolvent and repeat the 
same process. 

After the cleaning expansion, no inner tadpole terms nor counter-terms exist in the 
cleaned part in the dual graph. We shall check this nontrivial fact in section [6j However 

^The author is very grateful to V. Rivasseau for this suggestion. 
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the amplitude could still be divergent due to the possible arbitrarily large number of counter 
terms attached to the uncleaned part of the graph. Instead of further generating inner 
tadpoles we choose to re-sum them. After the resummation we would generate a divergent 
factor which we call the Nelson's factor. But this factor is not dangerous since the convergent 
factors we gained from the crossings and nesting lines would be enough to cancel it. We 
shall go back to this point in section [7} 

6 The renormalization 

6.1 The power counting theorem. 

Theorem 6.1. The amplitude of a graph of V vertices, B broken faces and genus g is 
proportional to 



where A is the ultraviolet cutoff J^, g is the genus and B is the number of broken faces. 

Remark that unlike the commutative model, the power counting for the vacuum 
Schwinger's functions is different from that for n-point Schwinger's functions due to the 
number of broken faces. 

We could prove the power counting theorem in at least three representations: the ma- 
trix base representation [5] [8], the direct space representation and [3l] and the parametric 
representation [35] . 

We shall not prove this theorem in this paper, as what we need essentially are the 
topological and combinatorial properties of the ribbon graphs. Since the LVE doesn't change 
the theory, it doesn't change the power counting either. 

6.2 The renormahzation 

Since the leaf tadpoles could be compensated by the counter terms before we use the LVE, 
the only possible divergences are the counter terms, and we use the inner tadpoles, the 
crossings and the nesting lines to compensate them. To be more precise, we use the inner 
tadpole to compensate the counter terms attached to the cleaned part and use the crossings 
and the nesting lines to compensate the uncleaned tadpoles. 

We have the following theorem for the renormalization, by which we mean the compen- 
sation between inner tadpoles and counter terms : 

Theorem 6.2. Each inner tadpole should be canceled by a counter-term. 

The compensation of the uncleaned tadpoles will be discussed in the next section. 

Proof. We consider an inner tadpole generated from an arbitrary resolvent R{a). We Taylor 
expand this resolvent as 




(70) 



R{a) 



1 



1 - iV2XCa + {iV2XCaf + Remainder, 



(71) 



1 + iV2\Ca 
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and the amplitude for an inner tadpole reads: 



Ttadpole = J dl^{a){y^jV2\CrnnCrnm)'^ = -2AT^. (72) 

mn 

The amplitude for a counter term attached to this resolvent reads: 



Tc = j du{cr)^^ X 2 X (iv^^^ Cmn(7nm)^rp(7ppTp = 2\T\ 



(73) 



So we have 

Ttadpole + T, = 0. (74) 

Remark that since the cancelation is exact, this theorem is true also for the sliced inner 
tadpoles and counter terms at any scale j. For proving this, we could just use the same 
method while replacing each propagator C by Ylij and each counter term Ta by Tj. 

So we have proved the remormalization theorem. We could also use the combinatorial 
renormalization (see the appendix) to prove the cancelation of the inner tadpoles with the 
counter terms. We don't discuss this point in detail in this paper as it is not essential. 
The interested reader could also look at [32] and |3T] for more details for the combinatorial 
properties for the LVE. □ 

7 Nelson's argument and the bound of the connected 
function 

In the last section we introduced the cleaning expansion, after which all inner tadpoles should 
be canceled by the counter-terms in the cleaned part of the dual graph. But there might still 
be arbitrarily many counter terms in the uncleaned part and they are divergent. Instead of 
canceling all of them, we resum them by using the inverse formula of the Gaussian integral 
and integration by parts. 

We write more explicitly, but loosely, the amplitude of the connected function after the 
LVE and the cleaning expansion as: 

^T= n [ Mw,a) f^w.TrfllKiMR^H^W^TA)], (75) 

where Rmn{<y) is the full resolvent that contains the pure propagator C, Tr needs to follow 
the cyclic order according to the real positions of the leaf terms K, the resolvents R and the 
counter terms Tj>^D. We have used the fact that all weakening factors for the counter terms 
equal to one, as they are leaves in the graph. There are only weakening factors for the a 
propagators that cross different regions in the dual graph. 

Now we consider the function G for which we haven't expanded the counter-term: 

0=1 du{a, w) W i:i[K^^Rmn{cj)]e''''^^"^\ (76) 

leT 
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We use the formula 



I dv{w,a)f{a)g{a) = e^l^^^-'-' -)^/(a)^?(a))U=o, (T7) 

where C{a,a',w) is the covariance that might depend on the weakening factor w or not. 
Hence 

G= [du{a,w)f2jr,[l~r{ n rrf^.Tr[irV)i?^"''^"^(-)]e^"^-^^]} 

N=0 ■ let ° 

OOOOOO „X rj 

Ari=0 Af2=0 Af3=0/Gr 

^^9(7'^ Ns\^2dada'^ ^''^=°- ^ ^ 

where we have written down expicitly the trace term in the exponential, and the term is 



defined by formula ( 20 ) . 



While the A^^i and N2 derivations generate connected terms, the last derivatives generate 
disconnected terms, see graph B in Figure [6] 

If we sum the the indices m for the counterterms directly, we would have: 

LAO O T-^,-.ATT„'T'. , t rt 



Af3=0 ^ N3=0 m=0 

= e-^^A., (79) 
And the resumed amplitude reads: 

A^ = l du{a, w) n Tr n[irL(^)^5.f (^)]e^^^^^'^e^^A, (80) 



.AT? 



We find that by resuming the counterterms we need to pay for a divergent factor e a, 
which is also called Nelson's factor. But as is shown in Formula (21), ~ A but not In A 



when A — )■ 00, so that this factor couldn't be compensated by the convergent factor e""*^^'^'^) 
that we have gained. 

For solving this problem we need to analyze more carefully the scale indices of the counter 
terms and other crossings before summing them |^ 

7.1 Improved Nelson's factor 



The reason for the bad factor A is due to that we sum all indices in ^m=o(^^'^^-^m)^ 



formula (79), which means we detach all counter terms in the uncleaned part of the dual 
graph even if they are convergent. For example, the tadpole of amplitude ln(A/m), where 
m is the scale of the border of the ribbon where the tadpole is attached, is not divergent if 



^The author is very grateful to V. Rivasseau for the discussion about this point. 
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m ~ cA, where c is a number that is not much smaller than 1, say 1/2. What's more, we 
need to bear in mind that each counterterm is attached to a pure propagator in the form of 
XCnmT^ and Cnm decays as M~^-^, where 2j ~ Inmax (m, n). Here we have ignore all the 
inessential factors. 

So instead of considering only T^, we need to take the whole object CnmT^, which we 
also call it the full counter term, into account. We first consider the case m > n, so we have 
m ~ M^-' where M > 1 is an arbitrary constant. Then we have 

C„^T^~-ln(A/m)<— . (81) 
m m 

So that the counterterm doesn't cause any divergence as long as m > In A and we could just 
bound these counter terms by 1 and need not to detach them from the connected graph. 
The counter terms become dangerous only when m < In A and we need to resumm them as 
introduced before. 



So that we only need to sum over the indices m for m < InA in formula (79), which 
reads: 



InA InA 



- 2A^ln2(A/m) < -2A^ln2A= -2A(lnA)^ (82) 

m=l m=l 

and the amplitude of the connected function after the resummation reads: 

A^= [ dHa, w) n Tr n[irL(c^)i?5.^ (^)]e^^^^^^e^(^°^)^ (83) 



whose divergent behavior is much better than in (80). 

If m < and < In A, we still need to resum the counter terms but the divergent factor 
should be smaller. If > InA, then the corresponding tadpole is convergent and we could 
just bound it by 1. 



This resummation process is shown in Figure 14 





1/2 3 

* exp[i (21) a T] exp[ IT ] 



Figure 14: A sketch of the resummation of the counter terms. Where the term T means InA. 
The counterterms with the scale index of outer border larger than In In A are not summed. 

The divergent factor e'^^^"^^^ is not dangerous as we could bound it by the following 
formula: 

as long as a is chosen properly, for example a = 1.4A. Here we have use the fact that fact 
that InA ~ and (ajL^)' ~ e^J^-i'^J™-. 
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So we have 



\AT\\T\=n < du{a,w)ll A"-i|Tr{||irLMII ll^if (^)ll le^-^'^^^lll 

X e^i'^'^e-'^^— ■ (ajLx)! ^ {KX)''-\ (85) 

Notice that the global trace operator Tr would also result in a cutoff index A, but this factor 
could be easily bounded by the exponential e""-'™"^ when a is properly chosen, and K is an 
arbitrary constant times the possible factors \/2 or \/2 + 1 comes from the bound of the 
resolvent R and the leaf K, respectively. 

This is the Nelson's argument in the formalism of Loop vertex expansion. 



8 Borel summability 



Theorem 8.1. The perturbation series of the connected function for 0| theory is Borel 
summable. 

Proof For the perturbation series Yln=o ^k^^ be Borel summable to the function G, we 
need to have 

G'(A) = ^afcA'= + i?("+^)(A), (86) 

n=0 

where R^'^^^\X) is the Taylor remainder. The analyticity domain C\ for G should be at least 
|A| < and ReA > [321 EH, which means 

-^<Argv^<^. (87) 

We rewrite the resolvent as 

R 



1 + ZV2AC1/VC1/2 

Since the matrix C^/^crC^/^ is Hermitian, its eigenvalues are real, so that there are no poles 
in the denominator. In the analytic domain of G we have 

II i? 11= I =^ 1^72, (89) 



and 

II is: 11 = 11 1 11^ 1 + v^. (90) 
However in the analytic domain C\ the linear counter term becomes: 



gTriV2A(TTA _ gTrj| V2A| cos ecrTAg-Trl V2A| sin ecrTA ^gj^^ 

where 9 = Arg-y/A. We could bound the first term in (91) by 1, but the second term would 
diverge for negative a. 
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We rewrite this term as: 



The term e'^^^^-^A could diverge at worst as e^^^^A for 9 = ±7r/4. But this is not dangerous, 
since we could still bound it with the convergent factor e~"-^™°^ ~ e~"^A that we gained from 
the crossings and nesting lines. 



We use simply the Taylor expansion with remainder for the connected function (86): 



Gmt=i = ^^l*=o+ f dt^^^^G^-^'\t\l (93) 



followed by explicit Wick contractions. We have for the reaminder 

||^n+i|| ^ |A|"+iK"(2n)!! < |A|"+i[K']"(n!), (94) 

where K and K' are positive numbers including the possible factors or \/2 + 1 from the 
bound of the resolvent Rmn and of the leaf Kmn respectively. Hence we have proved the 
Borel summability of the perturbation series. □ 

9 Appendix 

9.1 The order 1 case. 

There are only three terms for order 1 expansion and the corresponding graphs are shown 
in Figure [7} The amplitude for graph A reads: 

G\ = --[--(iV^A)'^] du{(7)^[{Crnn(TnpCrnpCrpn)] y- 3 

mnp 

= -In, (95) 

where the factor 3 comes from the number of contractions that can form graph A. The 
amplitude for graph B reads: 



Gb = ^(«v^)' j duia)[TrmT^ammTrnTnann\ = -XTl (96) 

Graph G is the constant counter-term with amplitude |T^. So we have 

G\ + G], + G'c = 0. (97) 
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B 






A. 



Figure 15: The renormalization of the divergent graphs at order 2. 



9.2 The order 2 case 



There are three kinds of graphs of order 2 which are shown in Figure 15 with nontrivial 
combinatorial factors. The amphtude of graphs A reads: 



G 



(2) 



T: 



A) 



1-{W2X) 



E 

mnlp 



1 



1 



m + p + ln+p + ll+p + 1 



X D 



(2) 



(9^ 



where 



E 

mnlp 



1 



1 



1 



m + p+ ln + p+ ll+p + 1' 



(99) 



(2) 

and D)^ is the combinatorial factor for graph A. For generating graph A we could either 
contract the half line 1 with 2 and 3 with 4, or 1 with 4 and 3 with 2. If the sigma fields 
are attached to the inner face like graph A2, these are all possible cases. If the a fields are 
attached to the outer faces, there are two additional cases for each kind of contraction, take 
Ai for example. So in total there are 8 possibilities and D\ = 8. Similarly we have = 9 



and D 



(2) 
C 



So the amphtude of graphs B reads: 



.(2) 
B 

3 



B 2! 



1 {iV2X) 

2 3~~ 



11 

^ ^ ^ ^n + p + ll +P+1 

nl 



X D 



(2) 
B 



(100) 



and the amplitude for graphs C reads: 



(2) 
C 



3! 



■3- 



1 Jn/2A)2^^^ 1 



X D 



(2) 
C 



(101) 
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So again we have 

Ga^ + + G*?^ = 0. (102) 
This checks exphcitly the renormahzation at low orders. 
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